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SUMMARY 


» 


The  governing  equations  for  laminated  beams  are  developed  by  variational 
methods.  The  beams  are  considered  to  consist  of  n  laminae  and  n-1 
bond  layers,  where  n  can  be  any  reasonable  number  greater  than  one. 

The  bond  thicknesses  are  assumed  to  be  small  in  comparison  with  those  of  the 
laminae.  The  effects  of  shear  strains  and  direct  strains  normal  to  the 
bending  axis  in  the  bond  layers  are  included.  Two  specific  problems  are 
solved:  bending  under  uniform  load  and  buckling  of  simply  supported 
beams.  Curves  are  presented  which  show  the  effect  of  bond  flexibility. 
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FOREWORD 


The  work  reported  herein  constitutes  a  portion  of  a  continuing  effort 
being  undertaken  at  Stanford  University  for  the  U.  S.  Army  Aviation 
Materiel  Laboratories  under  Contract  DAAJ02-68-C-0035  (Task  1F162204A17002) 
to  establish  accurate  theoretical  prediction  capability  for  the  static 
and  dynamic  behavior  of  aircraft  structural  components  utilizing  both 
conventional  and  unconventional  materials. 
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INTRODUCTION 


The  bending  of  multilayer  sandwich  plates  was  investigated  by  Liaw  and 

Little  in  Reference  1.  In  that  paper  the  authors  assume  the  facings  to 

carry  only  axial  loads  and  inplane  shear  (no  bending),  whereas  the  cores 

are  assumed  to  carry  only  transverse  shear  stress.  Stress  resultants  are 

then  formulated  which  reflect  the  contributions  of  the  cores  and  the 

facings  subject  to  the  above  assumptions.  Minimization  of  the  total 

complementary  energy,  in  a  procedure  similar  to  that  carried  out  by 
2 

Reissner  in  presenting  a  theory  for  the  bending  of  nonhomogeneous 
elastic  plates,  leads  to  the  governing  equations  for  a  multilayer  plate. 
Since  the  multilayer  effect  is  contained  in  the  stress  resultants,  the 
final  equations  for  the  stress  resultants  are  the  same  as  those  for  a 
single  sandwich  plate  except  for  the  definition  of  the  physical  constants. 

Inherent  in  the  Liaw  and  Little  analysis  is  the  assumption  of  the  same 

shear  angle  for  all  core  layers.  The  validity  of  this  assumption  was 

3 

investigated  by  Kao  and  Ross  ;  they  considered  a  beam  composed  of  three 
facings  and  two  core  layers.  By  allowing  the  structure  to  have  four 
degrees  of  freedom,  an  inplane  displacement  function  for  each  face  and  an 
out-of-plane  displacement  function  for  the  structure  as  a  whole,  the  two 
core  layers  are  then  capable  of  having  different  shear  angles.  Minimiza¬ 
tion  of  the  total  potential  energy  with  respect  to  each  degree  of 
freedom  leads  to  the  generation  of  the  governing  equations  for  the  system. 

In  both  References  1  and  3,  the  core  depths  must  be  considered  large  in 
comparison  with  the  thickness  of  the  facings.  Although,  in  Reference  3, 
the  bending  capability  of  the  facings  is  included,  the  strain-displacement 
relations  for  the  core  shears  are  not  valid  when  the  facings  are  thick 
since  they  do  not  include  the  effect  of  rotation  of  the  facings.  Refer¬ 
ence  1  provides  no  way  of  accurately  estimating  the  significant  bond 
stresses,  since  the  theory  deals  only  with  average  effects  across  the 
overall  plate  depth.  From  the  analysis  in  Reference  3,  an  estimate  of  the 
bond  shear  stress  is  possible;  however,  estimates  of  the  peel  strength  of 
the  bonds  are  not  possible  since  the  peel  strength  depends  on  the  direct 
stress  in  the  bond  layer  in  the  direction  normal  to  the  bending  axis. 


1 


Also,  In  Reference  3,  the  theory  end  its  application  are  limited  to  a 
beam  consisting  of  three  laminae. 

This  report  deals  with  laminated  beams  in  which  the  bond  layers  are 
small  in  comparison  with  the  thickness  of  the  laminae.  By  utilization 
of  variational  techniques  to  establish  the  governing  equations  and 
associated  boundary  conditions,  the  effect  of  the  bond  flexibility  on  the 
buckling  and  bending  of  laminated  beams  is  studied.  The  analysis  is 
performed  for  a  beam  consisting  of  n  laminae  and  n-1  bond  layers, 
where  n  can  be  any  reasonable  number  greater  than  one.  From  the  results 
of  the  bending  analysis,  the  shear  stresses  and  direct  stresses  normal  to 
the  bending  axis  in  the  bond  layers  can  be  calculated.  Such  information 
is  of  prime  importance  in  the  design  of  laminated  beams. 

The  laminated  beam  is  considered  to  have  2n  degrees  of  freedom,  an 
Inplane  and  an  out-of-plane  displacement  function  corresponding  to  each 
lamina  (see  Figure  1).  The  total  potential  energy  functional  is  then 
formulated  in  terms  of  these  2n  displacement  functions.  Coupling  of 
the  individual  laminae  is  achieved  by  describing  the  strain  energy  of  the 
bond  layers  in  terms  of  the  displacement  functions  of  the  laminae  they 
connect. 
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BASIC  THEORY 


STATEMENT  OF  PROBLEM  AND  BASIC  ASSUMPTIONS 

The  problem  considered  herein  is  that  of  analytically  formulating  the 
behavior  of  an  n-ply  laminated  beam  (see  Figure  1)  with  nonrigid  bond 
layers,  and  applying  this  formulation  to  the  problem  of  bend¬ 
ing  of  laminated  beams  under  a  uniformly  distributed  load  and  to  the 
problem  of  buckling  of  laminated  beams  under  axial  loading. 

The  development  is  predicated  on  the  following  assumptions: 

1.  For  the  individual  laminae,  all  of  the  assumptions  made  in  the 
engineering  theory  of  beams  are  assumed  to  apply. 

2.  The  laminae  all  have  the  same  elastic  and  geometric  properties. 

3.  The  bond  layers  all  have  the  same  elastic  and  geometric 
properties. 

4.  The  end  condition  for  the  structure  as  a  whole  is  that  of  simple 
support. 

5.  The  individual  laminae  at  the  ends  of  the  beam  are  not  restrained 
in  any  way  from  axial  movement. 

6.  The  bending  stresses  and  axial  stresses  parallel  to  the  bending 
axis  in  the  bonds  are  negligible. 

7.  The  stress  through  the  thickness  of  a  bond  layer  is  constant. 


STRAIN-DISPLACEMENT  RELATIONS  FOR  A  BOND  LAYER 


As  derived  in  Appendix  I,  the  shear  strain  in  a  bond  layer  (see 
Figure  2)  is  given  in  terms  of  the  axial  and  lateral  displacements  of 
the  i  and  i+1  laminae  as 
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STRAIN-DISPLACEMENT  RELATIONS  FOR  A  LAMINA 


In  accordance  with  elementary  beam  theory,  the  axial  strain  in  the  1th 
lamina  la 


€ 


X 


(2) 


STRAIN  ENERGY  IN  A  LAMINATED  BEAM 

The  strain  energy  stored  In  a  laminated  beam,  with. flexibility  effects 
Included,  is 

u  -  U1  +  u2  +  u3  (3) 

where  la  the  strain  energy  due  to  transverse  shearing  deformations 

of  the  bond  layers,  Is  the  average  strain  energy  In  the  bond  layers 

due  to  the  bonding  agent  compressibility  in  the  thickness  direction, 
and  is  the  strain  energy  associated  with  stretching  and  bending  of 
the  laminae.  The  quantities  and  are  derived  In  Appendix  I. 

Thus . 


VARIATIONAL  PRINCIPLE.  GOVERNING  EQUATIONS  AND  BOUNDARY  CONDITIONS 
Bending  Under  Distributed  Loading 

For  the  bending  problem  (see  Figure  3),  the  total  potential  energy 
stored  in  the  system  is 
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where  V  is  the  potential  of  the  distributed  load  acting  over  the  surface 
P 

of  the  i“  1  laminae. 

As  developed  in  Appendix ,11  the  extremes  of  U  +  with  respect  to  the 

2 n  degree*  of  freedom  u^  and  w^  (i  *  l,2,...,n)  are  the  2n  Euler 
equations 
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The  associated  boundary  conditions  obtained  from  the  variational  pro¬ 
cedure  carried  out  in  Appendix  II  are  the  3n  relations 
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Buckling  Under  Axial  Loading 

For  the  beam  buckling  problem  (see  Figure  4),  the  total  potential  energy 
stored  in  the  system  is 


where  Vp  is  the  potential  of  tne  axial  load  applied  over  the  end  «ur- 
faces  of  the  beam  laminae.  The  particular  relation  for  Vp  is  presented 
in  Appendix  III. 


As  developed  in  Appendix  III,  the  extreme  of  U  +  Vp  with  respect  to 
the  2n  degrees  of  freedom  u^  and  w^  (i  *  l,2,...n)  are  the  2n 
Euler  equations 
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The  Associated  boundary  conditions  obtained  from  the  variational  procedure 
carried  out  In  Appendix  111  are  the  3n  relations 
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METHOD  OF  SOLUTION 


For  the  two  problems  being  Investigated,  the  bending  of  a  laminated  beam 
under  lateral  loading  and  the  buckling  of  a  simply  supported  laminated 
beam  under  axial  loading,  the  governing  equations  and  boundary  conditions 
are  prescribed  in  equations  (6)  and  (7)  and  (9)  and  (10)  respectively. 

BENDING  UNDER  DISTRIBUTED  LOADING 

In  the  bending  problem,  the  complexity  of  the  equilibrium  equations  (6) 
precludes  straight forward  integration  and  solution  consistent  with  the 
boundary  conditions  of  simple  support.  Thus  w^  and  u^  are  assumed 
as 


w 


i 


i  =  1,2, . . .  ,n 
j  =  1,2 , ... ,® 


u 


i 


cos 


jftx 

L 


i  =  1,2 , . . .  ,n 
j  =  1,2,...,® 


(ID 


(12) 


and  used  in  conjunction  with  the  total  potential  energy,  equation  (5),  to 
effect  a  Rayliegh-Ritz  approximate  solution.  Both  u^  and  have 

been  selected  to  satisfy  the  geometric  boundary  conditions  for  simple 
support.  As  shown  in  Appendix IV,  the  application  of  the  minimum  total 
potential  energy  principle  leads  to  the  equilibrium  relations 


4  ■  f¥L  ■  G»b 


2t 


8L 


<*-W 


■ij 


E  bL  G  b  r,  .  * -2  ) 

_B_  +  _J_  IMeftlL  (1_&  J 

2t  t  8L  U  il;  i-1  j 


(continued) 
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jn  GBb(c?ft) 
4t 

r  j*  cBb(<^t) 


(1_8ii)*i-i  j  + 


i 
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jjt  GBb(eft) 


[  4t  J(1-5in>t+l  J  +  ["  2t  J(1‘5il)di-1  J 


4t 


GBbL 


1,2,. . .  ,n 

12  oo 

L  »  I  •  •  •  9 


<8i„-8nJ*ij 


Ste  iui!  +  Vi  .  .  8  , 

2  L  2t  1  il  in' 


U 


GBbL 

2t 


i  •  1,2, . . . ,n 
j  -  1,2,...,® 


In  matrix  form,  equations  (13)  and  (14)  become,  respectively. 


[Al]j  (ai)j  +  [A2]j  <Vj  *  <Vj 

[A3]j  (a^j  +  [A4]j  (di)j  -  0 


J  *  1,2, . . . / 


j  -  1,2,.../ 


(13) 


(14) 


(15) 


(16) 


The  matrices  [A^,  [A2 ] ^ ,  [A^,  and  [A^]^  are  given  in  Appendix  V. 


14 


Soludon  of  equations  (15)  and  (16)  for  (a^)  and  (d^)  gives 

"  [B^]1  (C^j  i  -  1,2,. ...n  (17) 

(dt)j  -  -[A^j1  •  [A3]j  [B^j1  (C1)J  j  -  1,2,. ...n  (18) 

where  [Bj]j  is  given  by 

[Bjlj  -  [Ax]  -  [Aj]  •  [A^j1  •  [A3]  j  *  1,2 . n  (19) 

The  application  of  equations  (17),  (18),  and  (19)  to  the  laterally 
loaded  beam  shown  in  Figure  4  leads  to  the  results  plotted  in  Figures  5, 
6,  and  7. 

BUCKLING  UNDER  AXIAL  LOADING 

In  the  beam  buckling  problem,  the  equilibrium  equations  (9)  and  the 
corresponding  boundary  condition  equations  (10)  are  all  satisfied  by 
the  displacement  functions 

wt  =  at  sin  ™  i  *  l,2,...,n  (20) 

-  d±  cos  ^  i  =  1,2,...  ,n  (21) 

Substitution  of  equations  (28)  and  (21)  into  the  equilibrium  equations 
(9)  leads  to  two  sets  of  n  algebraic  equations  as  shown  in  Appendix  V* 
In  matrix  form  these  equations  are 

[A5]  (at)  +  [Afi]  (dt)  =  0  (22) 

[A?]  (at)  +  [Ag  ]  (dt)  =  0  (23) 

The  matrices  [A,.],  [A  ],  [A^],  and  [Ag]  are  given  in  Appendix  V* 
Solving  equation  (23)  for  (d^)  and  substituting  into  equation  (22) 
gives 

[B2]  (at)  =  0  (24) 
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Figure  5.  Laminated  Beam  Geometries  Utilized  in  Bending  and  Buckling 
Solutions. 
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Figure  6.  Maximum  Deflection  of  Two-Fly  Beam  Under  Uniform  Loading. 


Figure  7.  Maximum  Deflection  of  Six-Ply  Beam  Under  Uniform  Loading. 


where 


[B2]  “  [A5]  -  [Afi ] [Ag ] ” 1  [A? ]  (25) 

For  nontrivial  solutions  of  (a^)  to  exist,  [E,,  ]  must  vanish.  Thus 
the  eigenvalues  of  [B^]  give  the  b<  ckling  values. 
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DISCUSSION  OF  RESULTS 


* 


The  equations  developed  in  this  report  for  the  bending  and  buckling  of 
laminated  beams  have  been  applied  to  two  specific  constructions,  one  a 
two-ply  and  the  other  a  six-ply  section  as  shown  in  Figure  5.  The  curves 
plotted  in  Figures  6,  7,  and  8  and  the  data  in  the  table  on  page  23  have 
been  derived  from  the  bending  solutions.  The  curves  plotted  in  Figure  9 
have  been  developed  from  the  buckling  solutions. 

Figures  6  and  7  show  the  effects  of  variations  of  the  bond  elastic 
properties  on  the  beam  deflections  for  the  two-ply  and  six-ply  beams, 
respectively.  The  deflection  values  have  been  nondimensionalized  with 
respect  to  the  deflection  obtained  when  the  bond  layers  are  infinitely 
rigid  in  shear  and  compression  normal  to  the  bending  axis.  As  shown  in 
Figures  6  and  7,  the  beam  deflections  are  significantly  affected  by  both 
the  shear  modulus  and  the  extensional  modulus  of  the  bond  layers.  For 
low  values  of  the  bond  shear  modulus,  the  laminae  tend  to  deflect  as 
individual  beams  and,  hence,  the  deflection  increases.  For  low  values 
of  the  bond  extensional  modulus,  the  bond  layers  tend  to  compress, 
causing  additional  lateral  deflection  of  the  beam.  There  appears  to  be 
little  quantitative  interaction  between  these  two  effects,  however,  as 
the  curves  in  Figures  6  and  7  are  essentially  parallel. 

The  shear  stress  in  the  bond  layers,  shown  in  Figure  8,  is  found  to 
depend  essentially  on  the  shear  modulus  of  the  bond  layers  and  to  be 
relatively  unaffected  by  variations  in  the  bond  extensional  modulus. 

The  shear  stresses  plotted  in  Figure  8  have  been  nondimensionalized 
with  respect  to  the  shear  stress  of  a  beam  with  infinitely  rigid  bond 
layers.  This  is,  of  course,  the  shear  stress  of  elementary  beam  theory, 
x  =  VQ/Ib.  As  shown  in  Figure  8,  the  actual  bond  shear  stress  does  not 
vary  significantly  from  the  shear  stress  of  elementary  beam  theory  for 
reasonable  values  of  the  bond  shear  modulus.  Furthermore,  it  is  useful 
to  note  that  the  elementary  beam  theory  shear  stress  is  conservative. 
Hence,  for  design  purposes,  the  elementary  beam  theory  shear  stress  may 
be  used. 
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Figure  9.  Buckling  Loads  for  Laminated  Beams. 


The  extensional  stresses  in  the  bond  layers  due  to  bending  are  given  In 
the  table  below  for  the  beaus  studied. 


MAXIMUM  DIUCT  SHKSS  II  BEAMS 

Eg/E 

o  /? 
max 

Two-Ply 

Six-Ply 

0.001 

0.001 

0.1 

0.00525 

0.00645 

0.00741 

0.00800 

0.00884 

0.00916 

*  ~ 

o,  max  bending  stress  in  beam  with  rigid  bonds 

The  extensional  stresses  in  the  bond  layers  are  found  to  be  relatively 
unaffected  by  variations  in  the  bond  shear  modulus.  Also,  as  shown  In 
the  table,  the  direct  Btress  is  not  greatly  affected  by  variations  in  the 
extensional  bond  modulus. 

The  buckling  loads  of  the  laminated  beams  studied  are  found  to  depend 
essentially  on  the  bond  shear  modulus,  with  variations  in  the  bond 
extensional  modulus  having  virtually  no  effect.  The  buckling  loads 
presented  in  Figure  9  have  been  nondimensionalized  with  respect  to  the 
buckling  load  of  a  beam  with  rigid  bond  layers.  As  shown  in  Figure  9, 
the  buckling  loads  are  significantly  affected  by  variations  in  the  bond 
shear  modulus;  in  fact,  for  beams  of  low  bond  shear  modulus,  the  laminae 
tend  to  buckle  individually.  It  is  of  Interest  to  note  that  although 
only  one  deflected  shape  is  assumed,  as  given  by  equations  (20)  and  (21), 
the  buckling  solution  yields  n  eigenvalues.  This  is  attributed  to  the 
fact  that  within  the  restriction  of  the  deflected  shape  implicit  in  (26) 
and  (21),  the  beam  can  still  develop  n-mode  ahapes.  That  is,  separation 
can  occur  at  one  interval  bond  layer  with  the  result  that  the  1  ami 
on  each  side  of  the  separation  move  in  opposite  directions.  Two  possible 
modes  for  a  four-ply  beam  are  shown  in  Figure  10.  However,  for  reasonable 
values  of  the  extensional  bond  modulus  (high  peel  strength),  the  eigen¬ 
values  corresponding  to  the  mode  shapes  shown  in  Figure  10  are  precluded. 
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CONC LUD ING  REMARKS 


The  analysis  presented  herein  is  the  first  step  toward  development  of  a 
general  theory  and  solution  procedure  for  laminated  plates  and  shells. 

The  distinguishing  feature  of  the  theory  is  that  the  transverse  shear  and 
extensional  flexibilities  of  the  bond  layer  between  two  laminae  are 
included.  The  extension  of  the  present  governing  equations  and  boundary 
conditions  to  plate  and  shell  elements  is  straightforward  with  the 
incorporation  of  the  appropriate  two-dimensional  strain-displacement 
relations  in  the  variational  principle.  The  quantitative  results  obtained 
for  the  bending  and  buckling  of  two-  and  pix-ply  beams  reflect  significant 
effects  of  bond  flexibility;  also,  they  indicate  that  redistribution  and 
amplification  of  the  bond  stresses  would  occur  under  conditions  wherein 
the  laminae  are  not  permitted  to  move  relative  to  one  another  at  the  beam 
boundaries.  This  certainly  is  the  case  to  be  expected  in  practice. 

Thus,  further  effort  should  be  expended  on  the  boundary  problem  as  well 
as  on  developing  both  qualitative  and  quantitative  information  for  laminated 
plates  and  shells.  In  qonclusion,  it  should  be  noted  that  two  of  the 
basic  assumptions  on  which  the  present  analysis  is  based  are  that  all  of 
the  laminae  and  all  of  the  bond  layers  possess  the  same  elastic  and 
geometric  properties.  However,  the  analysis  presented  herein  is  easily 
modified  for  beams  with  nonuniform  bond  layers  and/or  nonuniform  laminae. 
The  nonzero  elements  in  the  matrices  given  in  Appendixes  IV  and  V  can  be 
modified  to  accommodate  property  changes  between  the  laminae  and  the  bond 
layers . 
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APPENDIX  I 

STRAIN-DISPLACEMENT  RELATIONS  FOR  BOND  LAYERS  AND  STRAIN 
ENERGY  IN  BOND  LAYERS 


The  strain-displacement  relationships  for  bond  shear  and  extensional 
strains  normal  to  the  bending  axis  are  developed  in  this  appendix.  The 
consideration  of  two  adjacent  laminae  with  their  connecting  bond  layer 
and  the  distortion  of  a  typical  cross-section  a-a  (as  shown  in  Figure  2) 
leads  to  the  displacement  of  a  point  on  the  upper  edge  of  the  bond  layer 
due  to  bending  given  by 
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1 


c 

ui  ■  2  dx 


(26) 


At  the  lower  edge  of  the  bond, a  point  originally  on  section  a-a  is  dis¬ 
placed  a  distance 


.  c  dwi+l 
Ui+1  +  2  dx 


(27) 


Thus,  if  it  is  assumed  that  the  angle  of  distortion  of  the  bond  layer  is 
constant  across  the  thickness,  then  this  angle  of  distortion  is 


a  = 


V62 


c  %+l\ 

U1  Ui+1  2  \dx  +  dx  / 


(28) 


The  angle  a  includes  two  effects.  One  is  the  shear  distortion  of  the 
bond  layer  and  the  other  is  the  slope  of  the  bond  layer  due  to  bending. 
Under  the  assumption  that  bending  resistance  of  the  bond  layer  is  negli¬ 
gible,  this  latter  effect  is,  quantitatively, 


a 
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1  & .  fii±i\ 

2  \dx  dx  / 


(29) 


Thus,  the  angle  Q  can  be  written  as 


Use  of  equation  (28)  gives 


(30) 
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(31) 


u,  u  /dw.  dw.  \ 

=  -1-~ . -  —(1  +  ^)  (— i  +  — i±i 

t  2(  \dx  dx  / 


B 


For  thin  bond  layers  the  extensional  strain  normal  to  the  bending  axis  can 
be  assumed  constant.  Thus,  this  strain  component  is 


e  = 
z 


W  “W 

i+1  i 


The  strain  energy  in  the  bond  layers  is  then 


».-e|¥  ‘ 

i=l  ' 


}  |  .  c+t/d“i  , 

J  y  U1  ui+l  2  \dx  dx  / 
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(33) 
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APPENDIX  II 

GOVERNING  EQUATIONS  AND  BOUNDARY  CONDITIONS 
FOR  BENDING  UNDER  DISTRIBUTED  LOADING 

For  the  bending  problem  shown  in  Figure  3,  the  potential  of  the  applied 
loads  is 


V p  =  “  /  Pbwi 


The  strain  energy  associated  with  bending  and  stretching  of  the  laminae 
is 

v£hr /(S'*** /&)’-!  »■ 

i=l  v  0  0  'ax  '  ; 

From  equations  (33),  (34),  and  (35),  the  total  potential  functional  is 
n  /  -  L  ,  ,  v2  L  ,2  2 
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The  variation  with  respect  to  w^  gives 
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(continued) 
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Integration  by  parts  leads  to  the  set  of  n  coupled,  differential 
equations  and  2n  associated  boundary  conditions  giver,  by 
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The  variation  with  respect  to  u^  gives 
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Integration  by  parts  leads  to  the  set  of  n  coupled,  differential 
equations  and  the  associated  set  of  n  boundary  conditions  given  by 
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APPENDIX  III 


GOVERNING  EQUATIONS  AND  BOUNDARY  CONDITIONS  FOR 
BUCKLING  UNDER  AXIAL  LOADING 


For  Che  buckling  problem  shown  in  Figure  4,  the  potential  of  the  applied 
loads  is 

i  i  f  i(£)- 
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With  the  use  of  equations  (33)  in  Appendix  I  and  equation  (35)  in 
Appendix  II,  the  total  potential  functional  for  buckling  of  laminated 
beams  under  axial  loading  is  then 
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Integration  of  equation  (46)  by  parts  leads  to  the  set  of  n  coupled, 
differential  equations  and  the  2n  associated  boundary  conditions  given 
by 
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Integration  of  equation  (50)  by  parti  gives  the  let  of  n  coupled, 
differential  equations  and  n  associated  boundary  conditions  given  by 
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APPENDIX  IV 

SOLUTION  FOR  BENDING  OF  A  LAMINATED  BEAM  BY  A 
UNIFORMLY  DISTRIBUTED  LOAD 
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Varletlon  of  equation  (56)  with  respect  to  a^  leads  to 
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APPENDIX  V 

OLUTION  FOR  BUCKLING  OF  A  LAMINATED  BEAM  UNDER 
AXIAL  LOADING 

All  of  Che  governing  equations  end  boundary  conditions  developed  in 
Appendix  III  (equations  (47),  (48;,  (49),  (51)  and  (52))  are  satisfied  by 
the  displacement  functions 
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Substitution  of  the  functions  (67)  and  (68)  into  the  set  of  differential 
equations  (47)  as  given  in  Appendix  III  leads  to 
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The  set  of  equations  (69)  can  be  written  in  the  matrix  form 
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Substitution  of  the  displacement  functions  (67)  and  (68)  Into  the  set  of 
equations  (51)  In  Appendix  III  yields 
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